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With the start of the LHC, interest in electroweak scale models for the neutrino mass has grown. 
In this letter, we review two specific models that simultaneously explain neutrino masses and provide 
a suitable DM candidate. We discuss the implications of these models for various observations and 
experiments including the LHC, Lepton Flavor Violating (LFV) rare decays, direct and indirect 
dark matter searches and Kaon decay. 



I. INTRODUCTION 

o ■ 

Although the Standard Model (SM) of particle physics has been triumphant in explaining the observations in high 
energy accelerators, it has failed to accommodate the Lepton Flavor Violation (LFV) in the neutrino oscillation. 
Moreover, there is no suitable candidate for dark matter in the SM. Extensive literature exists on models explaining 
, each of these phenomena. For example, the conventional type I seesaw mechanism with very heavy right-handed 
■ neutrinos explains the small neutrino masses. However, this model does not provide a candidate for Dark Matter 
(DM). Moreover, since the scale of new physics in this scenario {i.e., the masses of the right-handed neutrinos) is 
much higher than the reach of the LHC or any other man-made accelerator in the foreseeable future, directly testing 
this scenario is a dream that may not ever come true. 

Testable models that can simultaneously explain neutrino masses and dark matter have been in the center of 
attention in recent years. To a great extent, this attention owes to the start of the LHC and a host of direct and 
indirect DM search experiments that are either in the process of data taking or planned to do so in the near future. 
To explain neutrino oscillation, the model should include sources for LFV. If the mass scale of the new particles is 
low, we in general expect rather large LFV decay rates, Br(// — » ej), Br(r — > ej) and Br(r — > ^7), that can be in 
principle observed at the experiments searching for these decays. With the MEG experiment taking data [l| and the 
superflavor project under consideration, this possibility is becoming more exciting. 

In this article, we review two specific models that can simultaneously explain neutrino masses and provide a DM 
candidate. Both models are based on a ^-symmetry under which all SM particles are even and all new particles are 
odd. As a result, the lightest new particle is stable against decay and provides a suitable DM candidate. In both 
these models, right-handed neutrinos are added to the SM; however, the Z2 symmetry prevents a Dirac mass term 
\ for neutrinos. Neutrinos acquire Majorana mass terms at one-loop level. 

One of these models, the so-called SLIM model, contains a low energy sector with masses less than O(10) MeV 
— ' [HOI- The low energy sector can show up at various observable phenomena like rare Kaon decay or supernovae. As 
t-H . we shall see, the upper bounds on the masses of these low energy sector particles and lower bounds on the couplings 
provide a way to test the SLIM model. The second model, the so-called AMEND model [J], does not contain such a 
• rH . low energy sector; however, it can lead to interesting phenomenology at direct searches for DM and at the LHC. 

The paper is organized as follows: In section [Hi we review the structure of the low energy SLIM scenario and its 
implication for low energy observations. In section ITOI we explain how the low energy effective SLIM scenario can be 
embedded within an electroweak symmetric model. In section IPVl we review the prediction of this model for various 
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observations such as DM annihilation into electron-positron or photon pairs, self-interaction of DM and LFV rare 
decays. In section IVl we discuss the possibility of discovering the model at the LHC. In section IVTI we first describe 
the structure of AMEND. We then review the neutrino mass matrix, LFV rare decay, DM annihilation and direct and 
indirect DM searches, electroweak precision test and signatures at colliders within the context of this model. Results 
are summarized in section fVIII 



II. THE REAL SLIM SCENARIO AND ITS IMPLICATIONS 

In this section, we review the so-called SLIM scenario which has been first introduced in Q ■ SLIM is the abbreviation 
of Scalar as Light as Mev which describes the characteristics of the DM candidate within this scenario. The scenario is 
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FIG. 1: One-loop contribution to the neutrino mass 

quite minimalistic and adds only a scalar, 5 and two (or more) right-handed Majorana neutrinos, Ni. A Z2 symmetry 
is imposed on the scenario under which S and iVj are odd but the SM particles are all even. As discussed in Q, the 
scenario can be realized both for real and complex 5. We will however concentrate on the real SLIM scenario. 6, 
which is called SLIM, is lighter than the other Z2 odd particles and as a result is stable and can play the role of the 
DM candidate. 

The scenario is based on the following Lagrangian: 



C 



(1) 



Of course, if we assume neither Ni nor S have electroweak interactions, the gi a coupling in Eq. ([T]) will not be 
symmetric under SU(2)xU(l). As a result, it can be only effective. That is, the scenario should be augmented to 
become symmetric under SU(2)xU(l) at higher energies. This step will be taken in the next section. In this section, 
we are concerned about the low energy implications of the scenario. Notice that the Majorana mass term for Ni 
violate lepton number and the coupling gi a is a source of LFV. As shown Q , this Lagrangian through the one- loop 
diagram shown in Fig. [TT] leads to the Majorana mass for active neutrinos 
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where A is the ultraviolet cutoff of the effective coupling which is the electroweak scale. Let us parameterize the 
neutrino mass matrix as follows 



m v = U ■ Diag[mi, rri2e 
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Combining (|2I3[) . we find 



g = Diag(Ai, X n ) ■ O ■ Diag( % /^T, \A^e 472 , y/m^e^U 7 
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where n is the number of Majorana neutrinos and O is an arbitrary n x 3 matrix that satisfies O t ■ O = Diag(l, 1, 1). 
Finally, 
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As is well-known, within thermal production scenario, the dark matter energy budget of the universe is almost inde- 
pendent of mass and is determined by the DM annihilation cross section. From this observation, the DM annihilation 
cross section is determined: (av) ~ 1CP 26 cm 3 sec~ 1 . Within the SLIM scenario the main annihilation modes are 
SS — > v a Vp,v a vp. Notice that these processes violate lepton number by two units. Indeed, these processes proceed 
through i-channel diagrams with a lepton number violating right-handed neutrino propagator proportional to : 
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Taking A — 200 GeV, 0.01 eV < m v < 1 eV and (<j(66 -> everything)?;) — 10 26 cm 3 sec _1 , from Eqs. (]2|6f> . we obtain 



0(1) MeV ~ m Nl ~ 10 MeV 



and 



3 x 10~ 4 ~ g la ~ 10" 3 



(7) 
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where Ni is the lightest right-handed neutrino whose propagator dominates the DM annihilation. Remember that 5 
is taken to be lighter than N so the DM candidate within this scenario is lighter than 10 MeV. This argument does 
not set a lower bound on mg; however from primordial nucleosynthesis a lower bound of O(l) MeV is derived 
That is within this scenario, the DM mass is in the MeV range. This is the rationale for naming 8 as Scalar as Light 
as Mev, SLIM. 

Although in this scenario DM pairs mainly annihilate to neutrino or antineutrino pairs, the energy of the produced 
neutrinos (E v ~ ms < 10 MeV) will be too low to be detectable at neutrino telescopes such as ICECUBE. In principle, 
as shown in Q, neutrinos produced by annihilation of 20-30 MeV DM candidates in dark halo can be observed by 
future large neutrino detectors such as LENA Q. However, for the mass range relevant for real SLIM, we do not 
expect a sizeable signal @. 

The most interesting feature of this scenario is that there is a lower bound on the coupling of the new sector to 
neutrino and an upper bound on their masses. This means by collecting enough data in the low energy observations 
that involve neutrinos, this model can be eventually tested. This feature and its phenomenological implications have 
been elaborated on in Q . Consider the decay A — > B + v where A and B can be detected but v appears as missing 
energy. If the mass difference between A and B is less than the sum of the masses of 8 and Ni, there will be another 
contribution to the missing energy signal: 

A -> BN8 

where both 8 and N, like v, escape detection. N eventually decay into is8. Thus, by studying the decay mode 
A — > B + missing energy, information on the parameters of the scenario can be deduced. Similar analysis has 
extensively been carried out (see Refs. fH-fHf) m the case of the Majoron couplings to neutrinos. In particular, consider 
the decay of K + to the charged leptons. By comparing T{K + — > e + + missing energy) /T(K + — > /i + + missing energy) 
with the SM prediction, the coupling gi a can be constrained: 

r( K+ -» e+ + missing energy) = T SM {K+ -» e+v e ) + ^ T(K+ -» e+N8) 

T(K+ -> M + + missing energy) ~ T SM {K+ -> n+v e ) + £ 2 T(K+ -> /i+JVtf) ( ' 

Tsm(K+ -» e+v e ) + EiT(X+ -> e+NjS) 
Tsm(K+ -> n+v e ) 

In the last line, we have taken 

L(A + -> e + N t 8) < L(A+ -> fjt+NiS) x 



L(A+ e+v) 

r(K+ -> fj,+v) 



Considering the fact that T(K + — ► e + z^)/L(A+ — > [i + v) ~ (m e /m^) 2 -C 1, this assumption is justified. The recent 
bounds on the ratio L(A + — > e + + missing energy)/L(A + — > /i + + missing energy) from KLOE [l3| yields 
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i 

where the sum runs over Ni lighter than K + . 

The spectrum of the charged lepton in the two-body decay K + — >• £ + + v will be of course different from that in the 
three body decay K + — > £ + + Ni + 8. Thus, by studying the spectrum of fj, + in K + — > ^ + missing energy, information 
can be derived on J2i \9ia\ 2 - The analysis based on 1973 LBL Bevatron data [l4| gives [|| 
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A more thorough investigation of K + — > n + Ni8 can be performed with the present KLOE data as well as with the 
upcoming NA62 results [15j . 

Another situation where the SLIM scenario can show up is the supernova explosion during which a neutrino gas 
of temperature of T ~ few 10 MeV is formed inside the supernova core. Since the SLIM particles are relatively light 
and are coupled to neutrinos, they can be produced at the supernova explosion. The produced 8 can interact with 
neutrinos in that environment with a cross section given by 0] 

Taking T ~ few x 10 MeV, we find that the mean free path of the SLIM particles is (crn^) -1 ~ 10 cm which is far 
shorter than the supernova core. As a result, the SLIM will be trapped inside the core. The energy transferred by 
diffuse out of the SLIM particles can be tolerated within the present uncertainties of supernova model • In case of 
supernova explosions in the future, the scenario might be tested by studying the neutrino energy spectra 0. 
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III. ULTRAVIOLET COMPLETION OF THE SLIM SCENARIO 



In the previous section, we introduced the low energy SLIM scenario based on the effective Lagrangian in Eq. ([T]) 
and briefly discussed its implications on the relevant low energy phenomena. The Lagrangian in Eq. ([1} has to be 
embedded within a SU(2)xU(l) symmetric model. In [2}, several ideas for the ultraviolet completion of the SLIM 
scenario have been suggested. In [3j a minimalistic model have been introduced that embeds the SLIM scenario. In 
this section, we review this model and in the next section we discuss its implications. 

The model presented in Q is quite minimalistic and is composed of (1) an electroweak singlet, 7/; (2) two Majorana 
right-handed neutrinos, iV,; and (3) an electroweak doublet with nonzero hypercharge, $ T = [<fi° <fi~] where <j>° = 
{4>\ + i(j) 2 )/V2 with real 4>\ and (f> 2 . As emphasized before, all these new particles are odd under the Z 2 symmetry. 
Imposing the Z 2 symmetry, the most general Z 2 even renormalizable Lagrangian involving only the scalars will be of 
form 
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Positivity of the potential at infinity puts constraints on the couplings [l6| , 

A1,A' 2 > 0, A 3 > -{K^) 1 ' 2 , A 3 > -(AA' 2 )V2 

and 

Ai -|a 2 | + a 4 >-(aa;; 
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Phases of A2 and m v ® are sources of CP-violation. For simplicity, we impose CP-symmetry on the Lagrangian in 
Eq. ([9]) which makes all the parameters in Eq. ((9]) real. 
Setting H T = (0 vh/V%), the mass terms will be of form 
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where 

0,2 

(11) 
(12) 

(13) 
(14) 

The parameters are taken in a range that neither of the scalars, except the SM Higgs, develops a vacuum expectation 
value. (f>~ , being a charged particle should be heavier than ~ 80 GcV to avoid the direct search bounds [17| ■ Notice 
that while (f>i mixes with r\ through the m ri § term, there is no such a mixing between <f> 2 and r\. Had we taken m, ( $ 
complex, the mixing term would be 

sR[m v $]v H <j) iri + ^s[m v q,]v H fori ■ 

However as we discussed above, we take the Lagrangian to be CP-symmetric so <fi 2 is a mass eigenstate itself. The 
other neutral mass eigenstates are 8\ and S 2 defined as follows: 

(15) 
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with 



tan 2a = 
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(16) 
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where in the last two equations we have taken (m r7 $Uiy) 2 /(m| i — m 2 ) 2 <C 1. In other words, the mixing, a is small 
and the interactions of lightest scalar Si with the W and Z-bosons are suppressed by sin a but 82 approximately 
corresponds to the real component of the electroweak doublet Direct searches [l8| restrict 82 and <j>2 to be heavier 
than ~ 90 GeV. On the other hand, 8\ can be light and play the role of the SLIM described in the previous section. 
To see this more clearly, let us add the couplings with fermions: 
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where L a is the lepton doublet of flavor a: = (i> La £ L ). We focus on the following range of parameters: 
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81, which is called SLIM, is the dark matter candidate. The main annihilation mode of DM is to neutrino (antineutrino) 
pair: 



(a(8i8i -> v La v L fi)v r ) = (a(8i8i -> v La v Lii )v r ) = 
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Considering that mg 1 < tuni , from this formula we expect the lightest right-handed neutrino to be one of the main 
contributors to the annihilation cross section so we find 
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Through a one-loop diagram active neutrinos acquire the following mass 
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This formula resembles the mass formula in Eq. @ with the difference that after UV completion, the UV cutoff 
has disappeared and instead the masses of the heavy particles, m^, 2 and m$ 2 show up in the formulas for the active 
neutrinos. For m 2 < m 2 N . <C m 2 , we find 



(m v ) a0 ~ 



32tt 



-m Ni 



sin 2 a{- 



log 



1 S1 



N, 



m 2 2 m%. 



■log- 



'jV, 



- 1) - A 



rr> 



Si 



2 — T 
ml 



(25) 



We can divide the parameter space to the following two separate regimes : (1) A2v|j / (mi ) ^> sin 2 a log(m 2 2 / rn? N ); (2) 
A2v|f/(m 2 2 ) ~ sin 2 alog(m 2 2 /m 2 v ) or Mvjj / (m 2 ^) <§C sin 2 alog(m 2 2 /m 2 vr ). In the first case, Eq. ([23| combined with 
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Eq. (|25l) implies mjVi "C 1 which is disfavored by big bang nucleosynthesis For A2i>ff/(m| ) ~ sin 2 a log(mj /m^), 
we find 

mjj <C mjVi ~ few MeV 

which is the same condition as we found for the low energy SLIM scenario in the previous section. 

Notice that within this model no upper bound on the masses of (f>~ , 8% or 2 is found. However, as we will see 
below, for relatively light $, the model is more natural. From Eq. ([25)1 and the perturbativity of gip, we find that 
sin a cannot be smaller than ~ 10~ 4 . For sin a <C 1 and m% 3> m 2 , we can write (see Eqs. (|16ll7p 

mg 1 ~ to 2 — sin 2 am% . (26) 

Thus, for m$ 3> 100 GeV, a fine tuned cancelation between the two terms in Eq. (|2"B)) is required to maintain ms 1 
below 10 MeV. In other words, based on naturalness of the model we expect <f>~ , 62 or 2 to be within the reach of 
the LHC. We shall discuss this point in section IVl 

Notice that this model has some features in common with the so-called inert model [l!| but with the difference 
that here we have an extra singlet scalar and the main annihilation mode of dark matter pair is into neutrinos. The 
contribution to the oblique parameters in our model is however similar to that in the inert model. Similarly to the 
inert model, within this model the SM Higgs can be as heavy as a few 100 GeV without violating bounds from the 
electroweak precision data. That is because the contribution from the new doublet to the oblique parameters can 
cancel the one from a SM Higgs. 



IV. IMPLICATIONS OF THE SLIM MODEL 



The impact of this model on the low energy phenomena such as the decay of light mesons and supernovae is similar 
to the low energy scenario discussed in section [TTI except that the coupling g^ has to be replaced by sina/v2- 
However, since this model also contains new heavy states, its phenomenology is richer. In particular, the heavy states 
can be produced at the LHC. We will discuss about this in more detail in section[V] Here, we discuss the impact that 
this model can have on other phenomena: (1) Annihilation modes of DM to an electron positron pair or photon pair; 
(2) Dark matter self-interaction; (3) Magnetic dipole moment of the muon; (4) LFV rare decay of charged lepton. As 
we shall see, bounds on rare decay already constrains a part of the parameter space. Let us discuss them one by one. 

A. Annihilation into electron positron pair and photon pair 

The annihilation to the e~e + pair is loop suppressed by a factor of e 4 /(167r 2 sin 2 9w) 4 0j- Because of this sup- 
pression, the rate of DM + DM — > e~e + is too low to account for the disputed 511 keV signal. Moreover the flux 
of radiation from the e~e + pair would be too low to be detectable. Of course, in this model the DM is too light to 
annihilate to fi~ fi + pair. 

At one-loop level a pair of 5\ can also annihilate into a photon pair with cross section 

ts r i e 8 sin 4 a m 2 - in -4i ( Ms t \ -4 3/ 

a(Si6i -t 77) ~ - 7— „.„ 7- t- ~ few x 10 ——7 sm 4 a cur/ sec . (27) 

v " ; 87r(167r 2 ) 2 cos 4 6W m^r V MeV / 

Because of the loop suppression, the flux of photons would be too small to be detectable at Fermi telescope (see, e.g., 
Fig 4 of [H). 

B. Self-interaction of Dark Matter 

The A' couplings in Eq. © can lead to the self-interaction of the DM pairs with the following cross section 

, . . . . . ,|A^| 2 sin 4 Q |A9| 2 cos 4 a IA0I 2 sin 2 acos 2 a, 
WWi hh)v) - Max L^U L_2J . 

The self-interaction of DM is constrained by merging galaxy clusters [2l| : a/rriDM ~ 1 cm 2 /g which translates into 

I A; 1 2 sin 4 a, I A 2 1 2 cos 4 a, | \' s \ 2 sin 2 a cos 2 a ~ 10" 4 . 
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C (9-2)„ 

Via coupling g^Ni^L^^) , muons receive a magnetic dipole moment at one- loop level: 

2 ^ 16tt 2 m 2 , V lJ ' 

where 



2g + 5^-1 gjogfr 
~ 12(^-1)3 - 2(^-1)4 ' (28) 



in which U = rn Ni /rn^ . For fj <C 1, 



s i^l ^ 5x 10 -i 2 E,l5 v l 2 f 100 GevV 



2 

which is two orders of magnitude below the present bound. 



D. LFV rare decay 

The coupling gi a Nit^ a (\r leads to the Lepton Flavor Violating rare decays, /i — > ej, r — > fj/y and t — » e7. From 
the formulas in [22j . we find 

777 

r(^^^ 7 ) = T7 fki ? l 2 , (29) 

1D7T 

where 

E* «em Q 
t 

where fj = rn 2 N Jm 2 ^ and K(ti) is defined in Eq. (|28|) . Within this model ii <C I however t2 can be either small or 
larger than f . In case that ta is also small we can write 

(4 
I^X) (30) 

Br(r 47) - 5 x 10" 5 | ]>> t5 *J 2 fl^X V . (31) 

s V / 

The latest bound [? ] on these branching ratios are 

Br(/i -> cy) < f .2 x lO" 11 (32) 
Br(r -> e 7 ) < 3.3 x 10~ 8 (33) 
Br(r -> M7) < 4.4 x I(T 8 . (34) 

These bounds already excludes large values of the couplings. However, the following range is consistent with bounds 
and is particularly interesting from phenomenological point of view as in this range the forthcoming LFV searches 
have a good chance of observing a signal: 

m^- ~ 100 GeV g ltJ ,,giT ~ few x 10~ 2 and g ie ~ few x 10 -3 . (35) 
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V. SIGNATURE OF THE SLIM MODEL AT THE LHC 



As we discussed in the previous section, naturalness implies that the new particles $T , 82 and <p2 have masses not 
much higher than 0(100 GeV). As a result, these particles are expected to be produced in pairs at the LHC. The 
produced 4>~ can decay via its coupling in Eq. (|19[) to charged leptons: 



M 2 ( 



m l- - m Nj 



T(4> -> LNi) = — 2 5 for m Ni < m*- - m\ B . (36) 

xbn rnr,_ 

Notice that the same couplings that determine the flavor structure of the neutrino mass matrix as well as branching 
ratios of the LFV rare decays of the charged leptons also determine these decays. This fact provides a means for the 
cross check of the model. In (24| , the feasibility of determining the coupling at the LHC has been studied considering 
the various sources of background and employing state-of-the-art techniques to enhance the signal to background 
ratio. In this section, we review the results of this analysis. The details of the analysis and the software that has been 
used to perform this analysis can be found in [24| . 

As we saw in section IIII1 at least the lightest N{ , which we call Ni has to be light which means the decay modes 
cf>~ — > N\e~ , Ni/i~ and N\t~ are all kinematically allowed. However, the second right-handed neutrino can be 
heavier. The following three situations are possible: 

• to tv 2 > TT!,^- . In this case, only (\>~ — > N\lg will be possible. 

• mjv 2 < and 777jv 2 ~ 777^-. In this case, both cj>~ — > Nxlg and (f>~ — > NzVj^ are kinematically possible 
and by studying the energy spectrum of the charged lepton, these two decay modes can in principle be distin- 
guished. However, considering that the contribution of N% to the neutrino mass cannot be larger than about 
~ \J Am 2 tm / Am 2 sun ~ 10 times the contribution of Nx, Eq. ([Ml) implies that 

for 777 N2 ~ 7770- , g 2 p < • 

This means for this situation, </> _ — > Nxl^ will dominate over 0~ —> N^ll. As a result, the signal for §~ — > 
lg + missing energy will be mostly composed of <f>~ — > Nxlp- 

• mjv 2 <C 7770-. In this case, the masses of both Nx and -/V2 can be neglected and the energy of lp in both cases 
will be approximately equal to in the rest frame of <j>~ . Thus, the signals for 4>~ — > Nxlp and <p~ — > N 2 lp 
cannot be distinguished. The signal for <f>~ — > la + missing energy is determined by the following sum 

\gip\ 2 + \92p\ 2 - 

In [24| only the case mjv 2 > 7770- is studied but from the above discussion, we conclude that the analysis in [24j also 
applies for the case itin-l *C nipf 2 < 7770- because 4>~ practically only decays to Nx just like the case m^ 2 > m^- . For 
the case ~ 777^, the analysis of [24| also applies but one has to replace \gxp\ 2 with \gip\ 2 + \g2p\ 2 - 

To make the analysis simpler it has been assumed in [24| that m0 2 — mj,-, ms 2 — m^- and |t77,5 2 — 77t0 2 | do not 
exceed 80 GeV to forbid two body decays 62, 4>2 — > W + cf>~~ or <p2 — > 82Z (or 82 — > 4>2%)- Moreover, <fi~ can have other 
decay modes such as <p~ W~8x, <p~ — > 8xt^v and (f>~ — > 8iW~j but these modes can be neglected for gx a ~ 0.01 
and sin a ~ 0.01. 

(23 | studies the pair production of (j) + 4>~ , (fi^fa and ^82 at the LHC and the subsequent decay of ^ to charged 
leptons. To perform the analysis, the two benchmark points with parameters shown in table Q] have been studied. At 
point A, BR^ ->■ e ± A^i) = and BR(0 ± fJ^Nx) ^ BR( ( /) ± r ± A^i) ~ 0.5. At point B, BR(</) ± e ± N 1 ) ~ 
BR( ( /) ± -> ^Nx) ~ BR(0 ± -> r ± iVi) ~ 1/3. The main focus in [H is on point A. 

The cross section of the 4> + 4>~ production at the 14 TeV run of the LHC for 7770± between 80 GeV to 130 GeV 
varies between 800 fb to 200 fb. At the benchmark point A, the subsequent decay of (jy^ will lead to four types of 
signal + missing energy, r + /i~ + missing energy, /i + r~ + missing energy and t + t~ + missing energy. The rr 

final state is contaminated by large hadronic backgrounds such as W + jets followed by the decay of W to jets or r. 
Light jets fake r even after applying the cuts so this mode has not been discussed in [24j . In [24| , the sum of signals 
t + fi~ + missing energy and /i + t~ + missing energy has been collectively studied. 

The main sources of background are W + W~ , ti, W+jets and Z+jets. Their cross sections at 14 TeV center of mass 
energy are shown in Table [TT1 Like the 4> + <fi~ signal, the W + W~ pair can lead to ltl~ + missing energy. Notice that 
the cross section of the W + W~ production is about two orders of magnitude higher than that of the + <?!>~ production. 
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Point A 


Point B 


m Nl (MeV) 


1 


1 


mjv 2 


> tn^- 


> ra^- 


a 


0.01 


0.01 


A 2 












m^ 2 (GeV) 




90 




90 


e 


tt/2 









f ° ^ 




1 0.01 \ 


9i<x 




0.03 




0.01 






V 0.03 / 




V -o.oi / 



TABLE I: Model parameters. 



Moreover, misidentification of some of the jets or other misidentifications can lead to mimicking the signal. To reduce 
the background and therefore enhance the signal significance (i.e., signal/ -^/background) , several cuts are suggested in 



24j . Using these cuts, explicit computation of the signal significance is carried out for benchmark point A at 14 TeV 
run of the LHC and for 30 fb _1 . The results are displayed for the r/x + missing energy and + missing energy signals 
respectively in tables UTT1 and ITVl As seen from table HVl for the values of m,*- lower than 130 GeV, the discovery 
can be made by 30 fb _1 of data. Considering that the background is almost the same, the signal significance of the 
e~e + + missing energy signal can be obtained by scaling that of pT /i + + missing energy by a factor of 

Br(0+ -> JV ie +)Br(<ft- -> Ni£)\m Point b _ 4 
Br(0+ -> A^+)Br(0- -> JVia*-)|a* Point A ~ 9 ' 

That is at m^i = 80 GeV and the benchmark point B, the signal significance of e + e~ + missing energy for 30 fb _1 of 
data will be 4.1a C.L. As indicated in [24]], this is a simplified estimation as there might be some difference between 
muon and electron reconstruction and selection efficiencies in the detector. A detailed study of these features needs 
a full simulation of the detector. 

A crucial question is whether g^ can be derived from the data. As discussed before, deriving the flavor structure of 
gip helps us to cross-check the model as the same couplings determine the neutrino mass matrix and LFV branching 
ratios. To derive gip, one should extract the signal number, N$ which is generally given by 

N S = N °^~ Nb (37 ) 

where 65 is the selection efficiency of the signal, V b s . is the observed number of events and Ng is the contamination 
due to the background. Nb is calculated by simulation using input such as Parton Distribution Functions (PDF) or 
total luminosity. The uncertainty in these inputs induce an uncertainty in extracting N$ and therefore the couplings. 
The main source of uncertainty is the uncertainty in PDFs which at present is about 10 % uncertainty. These 
uncertainties induce an uncertainty of about 60 % in extracting Ns- In [23], it was shown that ANs/Ns due to 
these uncertainties does not improve by increasing the luminosity. However, increasing the center of mass energy will 
enhance signal to background ratio and therefore improve ANs/Ns- Of course, if by using the data of LHC or some 
other machine, the uncertainty in PDFs are reduced, the precision of extracting gi a can be improved. 

Other modes that have been discussed in [24| are pp — > </> ± ^2 and pp — > 4> ± S 2 . As discussed earlier in this section, 
because of the simplifying assumptions on the mass spectrum of the components of <&, decay modes such as (j>~ — > W~5\ 
are forbidden so 62 and (j) 2 can have only invisible decay modes 62, 4>2 — > N\U. As a result, the signal will be composed 
of a charged lepton from the decay of (j)^ plus missing energy which is composed of Ni from ^ decay and the decay 
products of 4> 2 or 62. At benchmark point A and for m^± = 80 GeV, the significance of the /1 + missing energy and 
t + missing energy signals can reach as high as 9er and 4.6<r at 14 TeV energy and 30 fb _1 integrated luminosity [24j . 
Again under simplifying assumption that the selection efficiency of detecting muon and electron is not much different 
the significance of the e + missing energy signal at the B point is equal to that of /x + missing energy signal at the A 
point rescaled by a factor of Br(0+ -> Vie + )| A t point B/Br(^+ -» Vi^+)| A t point A ^ 2/3. 

The above results are for the 14 TeV run of the LHC. In [24], an estimation of the signal significance for the 7 TeV 
is made by rescaling the cross sections of both background and signal to their values at 7 TeV run. The results for 30 
fb _1 are displayed at table IVl 
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Process 


W + W~ 


tt 


W+jets 


Z+jets 


Cross Section 


115.5±0.4 pb 


878.7±0.5 pb 


187.1±0.1 nb 


258.9±0.7 nb 



TABLE II: Background cross sections 





80 GeV 


90 GeV 


110 GeV 


130 GeV 


Signal significance 


2.8 


2.2 


1.4 


1 



TABLE III: Signal significance in the TfiE™ las final state for different m,<±^ hypotheses at 30 fb 1 



VI. AMEND 

In this section, we review the AMEND model which was introduced in [4j . AMEND stands for A Model Explaining 
Neutrino masses and Dark matter. Like the SLIM model, there is a Z2 symmetry that protects DM from decay. The 
Z2 symmetry also forbids Dirac mass at one-loop level. The Z2 symmetry in this model is the remnant of a global 
Xl(l)x symmetry which is explicitly broken by small parameters. In the limit of exact U(l)x, neutrino masses vanish. 
Neutrino masses are suppressed both by a loop factor and the small U(l)x breaking term (i.e., 't Hooft criterion). The 
particle content of the model includes two fermionic doublets R and R' with opposite hypercharges, an electroweak 
triplet A and a complex singlet, cf>. These particles can in principle be produced at colliders. In particular, one of 
the components of A is doubly charged and can lead to almost background free signal of same sign charged lepton 
pair plus missing energy via A ++ — > l^ltS^. The LFV couplings of these particles lead to LFV rare decays such as 
fi — >• &y, t — y e^f and r — > fxj at loop level. 

Within this model, the DM production in the early universe is thermal. In [3] , various possible annihilation modes 
have been discussed. It was found that the dominant annihilation mode was the DM annihilation through s-channel 
Higgs exchange which can account for the observed density of DM. DM in this model can be counted as a Weakly 
Interacting Massive Particle (WIMP) and can show up in direct DM search experiment based on measuring recoil 
energy from scattering DM off nuclei in a background free environment. Various experiments are designed for this 
purpose. Their results are contradictory. On one hand, experiments such as XENON do not find any signal and 
on the other hand, the DAMA experiment [25[ reports a positive signal at more than 8 <r C.L. To reconcile these 
conflicting results, several attempts have been made. Among these solutions, inelastic DM solution [27| and light 
DM (< 10 GeV) [iH have received more attention. In null results interpreted as an upper bound on the cross- 
section as well as the two solutions accommodating the positive signal from direct DM searches have been studied and 
shown that they can be embedded within the AMEND model by going to the proper regions of the parameter space. 
However, more recent data from XENON100 [29j and a re-analysis of XENON10 data [26| respectively disfavor the 
inelastic DM and light DM solutions. We therefore focus on the constraint from DM searches. 

In this section, we first describe AMEND and then discuss its implications for various observations. 



A. Description of AMEND 



The particle content of this model is listed in table lVIl As seen from this table, the new particles are the following. 

• A complex scalar field which can be decomposed in terms of real fields as <j> = (4>i + i<f>i)/ V2; 

• A scalar triplet with the following components 



A = 



Al A++ 



■J2 
A 



A + 
' V2 



(38) 





80 GeV 


90 GeV 


110 GeV 


130 GeV 


Signal significance 


9.2 


8.4 


6.6 


4.9 



TABLE IV: Signal significance in fi/j,E™ l3a final state for different mr^-i hypotheses 
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Channel 


Mass Point 


Signal significance 




m w±) = 80 GeV 


1.6 


(j> + <f>- -> t^ET ss 


m w±) = 90 GeV 


1.2 




m (0±) = 110 GeV 


0.7 




m {4>±) = 130 GeV 


0.5 




m (0 ±) = 80 GeV 


6.4 


<j>+<j>- -> n)iE^ sa 


m w±) = 90 GeV 


5.7 




m (0±) = 110 GeV 


4.2 




m w±) = 130 GeV 


3.0 


-> t£^ ss 


m (0±) = 80 GeV 


2.6 


4> 4>2 ^ (j,e™ tas 


m (0±) = 80 GeV 


5.0 



TABLE V: Signal significance in different final states for the 7 TeV run, provided that 30 fb 1 data is collected at this energy 
before any switch to higher machine energies. 



particle 


SU(3) C SU(2) L 


U(1) Y 




Ql 


3 


2 


1/6 




Ur 


3 


1 


2/3 




d R 


3 


1 


-1/3 


fermion 


II 


1 


2 


-1/2 




&r 


1 


1 


-1 




R = Rr 


1 


2 


-1/2 




R' = R' R 


1 


2 


1/2 




H 


1 


2 


1/2 




A 


1 


3 


1 


scalar 


<t> 


1 


1 








TABLE VI: Particle content and gauge quantum numbers. 



We can write A = (Ai + iA2)/V2 where Ai and A2 are real scalar fields. 

• Two Weyl fermion SU2 L doublets, R T = {v R E R ) and (R') T = v' R ). 

A symmetry, called G symmetry, is defined under which each of the new particles are charged under a separate U(l). 
The G symmetry is defined as follows. 

G = U(l) K xU(l) xU(l) A xU(1) ( (39) 

where U(l)^ and U(1)a are symmetries under which only <j> and A are respectively charged and U(1)_r is the symmetry 
under which R and R' have opposite quantum numbers. U(l)^ is the familiar lepton number U(l) symmetry associated 
with lepton number. The model is constructed such that the main part of its Lagrangian preserves the G symmetry. 
In addition to the kinetic and the gauge interaction terms, the most general G-preserving Lagrangian is composed of 
the following scalar potential 



V = - n%H*H + n 2 A Tr (At A) + ^<f> 

+ \(H^Hf + ^<j>) 2 + ^ ( TrAtA) 2 + ^ Tr(At[At, A] A) (40) 
+ XhaiH^H Tr (At A) + X HA2 H^[A^ , A]H + Tr (At A) + \ H4> tf <j>H* H , 

and the fermionic part which is a Dirac mass term for the fermionic doublet 

-C R = m RR {R' c )^ -R + h.c. , (41) 

where (R' C ) T — {v' R — (E R ) C ). m RR should be heavier than ~ 100 GeV to avoid the bounds from direct searches. 
As a reference point, we shall take m RR — 300 GeV. At a very high energy scale, the G symmetry breaks to a smaller 
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U(l)x symmetry under which the SM particles are all neutral and the quantum numbers of the new particles are as 
follows. 

D % , 1 TV U ili x 1 . U(l\x . 1 , , U{l)x . 

R => +1 , R => —1 , A => +1 and <p => —1. 

Notice that U(l)x is anomaly-free and can in principle be fixed. The terms that break G to its U(l)x subgroup are 
the following 



V. 



HAd 



--X 



HA(f> 



H T \a 2 I^H(^ + h.c. 



Vi L 4> =g a <j> ] R 1 '4a + h.c. 



(42a) 
(42b) 



At a lower energy scale, U(l)x breaks to a Zi symmetry under which the SM particles are even but the new particles 
are all odd. After U(l)x —> Z%, the Lagrangian includes the following terms for the scalars 

Vscalar = X H A<t>H T ia 2 A^H<j) + p?^ 2 + X<p ^ + A* 2 ^V + Xh^Hc/) 2 + X A ^tA^ A(j) 2 + h.c. . (43) 

and the following terms for the fermions 

- C eL<t> = g a cf>R ] £ La + h.c. and - C £lA = (sa)^' 1 ■ A ■ £ La + h.c. (44) 

The pattern of symmetry breaking implies g ^> g,gA and Xha<p 3> ^HA4>- 

B. The scalar sector 



Within this model, only the SM Higgs obtains a vacuum expectation value and 

(<M = (<h) = (Ai> = (A 2 > = . 



(45) 



After electroweak symmetry breaking, <j> and A° mix with each other. For the CP-symmetric case, CP even scalars 
4>i and Ai mix only with each other and the CP-odd scalars <f>2 and A2 mix among each other. The neutral scalar 
mass eigenstates, 8\, <5 2 , ^3 and 64 can be written in terms of the components of <f) and A as follows 



(46) 



where |tan2ai| ~ |tan2a2| — 2m| A /(m^ — to?). In the following the masses of Si are denoted by Mj. The formula 
for Mi can be found in The difference ||ai| — 1 0:2 1 1 as well as the mass splittings \M% — Mi| and \M& — M%\ are 
suppressed by the U(l)x-breaking terms. We take 8\ to be the lightest new particle and therefore a DM candidate. 
As discussed in Q, the CP-odd scalar S2 could also play the role of DM. 
The coupling of 61 to the Z-boson is of form 



(sA 




cosai 





sinai 


\ 




s 2 







COS OC2 





sin a.2 


4>2 


s 3 




— sin o.\ 





cosai 





Ax 




\ 





— sin ot2 





COS«2 / 


va 2 ; 



i5su(2) smai sma 2 



cos 9w 

where 5su(2) is the SM weak gauge coupling and 9w is the Weinberg angle. If Mi + M 2 < mz, 



Gf sin ai sin 012 
6y/2n ' 



(47) 



(48) 



82 will eventually decay to S\ and a neutrino pair so this decay mode will count as an extra contribution to invisible 
decay mode of the Z boson. For Mi + M2 < mz, the upper bound on the extra invisible decay modes of the Z boson 
23] implies 



sinai sina2 < 0.07 
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C. Neutrino masses and LFV rare decays 



Within this model, there are one-loop contributions to the neutrino mass matrix of form 

K)a? = [ffa(<?A)/3 + 9p{$A)a]V + [SoiGm.)/? + 5/3 (<?a) q ]?? , 



(49) 



m RR 
64tt 2 



Mf 



rn 



RR 



Mf 



In 



"RR 



Mf 



Ml 



rn 



RR 



Ml 



In 



lb RR 

Ml 



sin2ai - [{ax, Ml, 



Ml) 



a 2 



Ml Mi)] , (50a) 



m RR 
64tt 2 



Mf 



"RR. 



Mf 



In 



Mf 



Ml 



m 



RR 



Ml 



In. 



11 RR. 

Ml 



sin2 ai + [(«!, M 2 , M|) -> (a 2 , Af|, M 4 2 )] . (50b) 



The parameters denoted by tilde are all suppressed by U(l)x breaking terms. Notice that for g A = 0, the neutrino 
masses vanish. This is expected as for g A = 0, by assigning lepton number equal to +1 and —1 respectively to R and 
R', lepton number will be conserved so the neutrinos cannot have a Majorana mass term. It is straightforward to 
confirm that regardless of the flavor structure of the couplings, the determinant of m„ vanishes which means one of 
the neutrino mass eigenvalues is zero and the neutrino mass scheme is hierarchical. This structure is due to the fact 
that only two right-handed neutrinos are incorporated within this model. In order to make the neutrino mass scheme 
non-hierarchical (i.e., Det[m„] ^ 0), another pair of R and R' should be added. 
For m RR = 300 GeV and m v = 0.05 eV, it has been found [3| that 



99A 



3.4 x 10" 



70 GeV 50 MeV 



TUrr 



0.1 



0.05 eV Mi 



300 GeV I sin I 



L RR. 



■In- 



L RR 



rn 



RR 

for 



+ 1 - In ■ 



IL RR 

Ml 



2m^m^ A /m A ~ 2Mi<5| sin a x | > m% A , (51a) 



gg A ^ 3.3 x 10" 



300 GeV 1 GeV' 



0.05 eV m RR 



™\a 



f toa y 

V500 GeV/ 



l RR. 



in. 



A 

for 



In- 



2m\m\ A 



^,2 



2Mi5\ sin ai| < rh 



<A A 



(51b) 



gg A ~ 1.3 x 10" 



300 GeV 0.1 



"RR 



0.05 eV m RR jsinai 



m A f ln m RR 



(52) 



The gp, gp and (<Ja)p couplings will lead to LFV rare decays such as l a — > £97 4j. Since (<?a)/3 , gp <C gp, the 
dominant contribution is from the g coupling: 



Br(^i ->■ ei) « 2.5 x 10~ 9 
Br(r l al ) w 4.5 x 10 



/300 GeV x 


4 


5^ 9e 


V m RR / 


) 


0.10.1 



and 
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/300 GeV\ 4 


9*r 9a 


V m RR J 


0.10.1 



(53a) 
(53b) 



Let us now discuss the constraints on the parameters from bounds l a — > Ipj [? ]. The bounds on Br(r — > ej) and 
Br(r — > fi-f) allow even values of m RR as small as 100 GeV and <7 MiT as large as 0.2. For g e , g^ ~ 0.1, the bound on 
Br(/i — > e7) requires relatively large values of m RR , m RR > 1.1 TeV. However, for g^ ~ 0.02 and g e ~ 0.01, as 
small as 100 GeV can still be accommodated. An alternative solution is g e <C g^ or g e g^. In the case g e <C <7 M , 
the ggA'H contribution dominates (m v ) ea ; i.e., (m v ) ea = [g e (gA)a + <7a(<?A)e]?7- Similarly for the case g e 3> g^, 
(m v )^ a = [g^{g A ) a + 9a(gA) l ,}v- 



D. DM annihilation and searches for DM 

As discussed before, within this model the DM production is thermal. To obtain the observed amount of the DM 
density, the annihilation cross section should be 

(cr(5i<5i ->• anything)^) ~3x 10~ 26 cm 3 /sec , (54) 
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where v is the relative velocity. In fact, since the mass splitting between 5% and 82 might be small, the effect of 82 at 
the freeze-out epoch has to be taken into account. 

In [4|, different possible annihilation modes were investigated. Depending on the regions of the parameter space, 
different annihilation modes dominate. These modes are list as follows: (1) Higgs mediated decay into //: <5i<5i — > 
h* -> //; (2) Higgs mediated decay into an on-shell W plus and an off-shell W: 8\8i -> WW* -> Wff'; (3) 
Annihilation into WW or WW* via gauge interaction; (4) Annihilation into a Higgs pair. Of course annihilation into 
W or H pair can bepossible only for heavy Mi. The annihilation of these particles in the sun center can lead to a 
hard neutrino flux [30] which is disfavored by the present indirect dark matter searches. Annihilation modes to rf 
and cc are also disfavored by indirect DM searches at the neutrino telescopes. We shall focus in the range for which 
m& < Mi < 70 GeV < m^. In this range, the dominant annihilation mode is S-^Si — > h* — > bb so the emerging 
neutrino flux from DM will be rather soft and can be tolerated within the present bounds. 

The coupling of 82 and 8\ to the SM Higgs is given by 



XLVnhSf = -y ((Ajjai - X H & 2 ) sin 2 a% + \h<p cos 2 a\ - 2A ff A0sinai cosai) hS { 

(55) 



Mf — [i^ cos^ a,\ — ^ sin a\ 



hS: 



2 



where i = 1,2 and the sub-dominant U(l)x violating terms are neglected. Notice that the couplings of U(l)x 
preserving part of the Lagrangian can be made real by redefining the fields. In particular, Xha<p can be made real so 
there will be no coupling of type hS\S2- The coupling of form h8\82 appears when both CP and U(l)x symmetries 
are broken by e.g., 3[AijA<^]- 
The Xl coupling leads to 

- v |A L | 2 m 2 (M 2 -m 2 ) 3 / 2 

where nif is the fermion mass for the kinematically accessible channels and N c — 3 (1) for quarks (leptons). In the 
limit (M2 — Mi)/(2Mi) <C 1, we have to take into account the annihilation of 82S2 in the calculation of the DM 
abundance as 

(a(5 2 S2 -> ff)Hv) ~ (<t(6i6i -> ff) H v) . 
Eventually 82 decays via coupling to a virtual Z boson exchange Q : 

, -s 1 . (M 2 - Mi A 5 /sinai\ 4 _ x 

T(8 2 -> ~ 15 sec . 

v 2 1 ; V 50MeV y V 0.1 J 

We therefore generally expect the decay to take place before Big Bang Nucleosynthesis (BBN) so it should not affect 
the BBN predictions. 

The coupling in Eq. ([53]) leads to the interaction of DM with nuclei via i-channel Higgs boson exchange so it can be 
constrained by direct DM searches at underground experiments sensitive to the recoil energy of the nuclei interacting 



with DM particles. The cross section of the DM with nuclei can be written as [31 1 



I A, I 2 « /2 ^ 6 5 x 1Q _ 45 (Xr_Y f 65GeV\ 2 (l20GeV\ 4 U . ^ ^ 



p -n Mlm\ J V0-047 V M i J \ m h J \0.2 

where ^s in is the reduced mass of the dark matter- neutron system, m n is the nucleon mass and / is the nuclear 
matrix element parameter which can vary within the present uncertainties, 0.14 < / < 0.66 [3l| . The same coupling 
also determines annihilation of DM so we can write 



/ M 2 m m p( 4M i - m h) (<r(Si5i -> h* -> SM final states) v) 
irM-Lmfrvjj 4T(h -t SM final states) | mfe _>2Mi 



(58) 



The strong bound on [32j already constrains a considerable part of the parameter space. We generally expect a signal 
at future DM search experiments unless Mi — > m/j/2. The case of Mi — > m^/2 seems to be unnatural. If the future 
searches for DM do not find a signal and no neutral stable scalar with mass equal to m^/2 is found by colliders, this 
model will be disfavored and will eventually be ruled out. 
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E. Electroweak precision tests 



The new particles added to SM participate in the electroweak interactions so they can lead to corrections to 
the electroweak precision parameters. The contributions to the oblique parameters S, T, W, Y (33j were explicitly 
calculated in Q- 

Since R and R' have equal masses, their contributions to S and T parameter exactly cancel each other. Their 
contribution to W and Y will be non-zero but tiny: 

120n 2 m 2 RR 120tt 2 m\ R 

The contribution of A to the oblique parameters can be written as 

25 0qtt/ o % 777,2 7<7qtUT\ TJ1.?«, 7Q? T/1 \ 777?., 



a _ -9su(2) a ffsp(2) mj 2 [ ffsu(2) ' 9 v{1) m z w 

24^ 2 ? ' 576tt 2 m 2 / ' 720tt 2 to| ' 

where the relation 2 to* + = to 2 ^ + has been used and the results have been expanded in 



+ -"i =AgAa 4. (60 ) 

The contribution of </> to the electroweak precision parameters is suppressed by a factor of | sin a\ sin «2 1 relative to 
that of A so it can be neglected. Within this model, Higgs can be heavier than in the SM model because the new 
contributions can cancel out the contributions from the Higgs to the oblique parameters and the upper bounds from 
electroweak precision tests on the Higgs mass can be relaxed. Without cancelation {i.e., for a light Higgs mass), the 
T parameter constrains £ < 0.1 which translates into a bound on the splitting of the components of the triplet. This 
results in a mild bound on XhA2, e.g., for toa — 500 GcV, the bound is XhA2 ^ 0.5. 



F. Signatures at colliders 

We expect a rich phenomenology for the LHC within this model. For M\,M2 < mh/2, the SM Higgs boson can 
decay into a pair of Si or a pair of 62- These new decay modes can dominate over the decay to the 66 pair when 
Xl ~ rrib/vH — 0.02. Decay mode to Si pair will appear as missing energy. If AI2 — Mi is larger than 2m e , 62 can 
also decay into S\e~e + which appears as a distinct displaced vertex. In order for the decay to take place within the 
detector, the following condition is necessary: dTs 2 /2"f > v where v is the velocity of 62 and 7 = (1 — u 2 ) -1 / 2 . d 
characterizes the size of the detector. As shown in this condition requires \M.2 — M\\ > 500 MeV. Remember 
that M2 — M\ is suppressed by the U(l)x violating parameters. For smaller values of the splitting the decay mode 
H — > S262 will also appear as missing energy signal regardless of if the 62 ~ > 8\e~e + mode is kinematically accessible 
or not. 

If the new particles are not too heavy, the charged particles A ++ , A + , E R and E R as well as the neutral particles 
#3, <$4, Vr and v' R can be produced through electroweak interactions. They will eventually decay into the SM particles 
plus Si or ($2- In particular A ++ can decay to a pair of same-sign charged leptons: 

A++ -> l+l+Si, 2 ■ 

The background from the SM to the same sign charged lepton pair signal is not very high which makes the discovery of 
A ++ easier. T(A ++ — » £%£g Sip) is proportional to \{g&)agp + {g&)pga\ 2 ■ The decays of the charged components of R 
and R' are given by the g a couplings: Br(E R — > £^81,2) cx \g a \ 2 . In principle, g a and (§a)<x can be directly extracted by 
studying the flavor pattern of the decay modes of A ++ and Er. Notice that the coupling determining Br(E R —> £~Si t 2) 
also determines Br(Z~ — > ^7). Moreover the combination {gA) a gp + (5a)/35q determining T(A ++ — > £££p5i,2) is 
exactly the combination appearing in the neutrino mass matrix in Eq. (|49j) . 



VII. CONCLUDING REMARKS 

In this letter, we have first reviewed the SLIM scenario and the model embedding it p|. We have then reviewed 
AMEND, which stands for A Model Explaining Neutrino masses and Dark matter [4]. Both these models are con- 
structed to explain the tiny neutrino mass and provide us with a DM candidate. We have reviewed the implications 
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of these models in various experiments and observations. We have updated the results in [2|-|4(, taking into account 
the most recent data release from experiments such as direct DM searches XENON10 [26| and XENON100 [29j . 

In both models, there is a Z2 symmetry that makes DM stable and forbids a Dirac mass term for neutrinos. The 
neutrinos acquire Majorana masses at one loop level. In the framework of the both models, we generally expect the 
value of Br(/j, — > ej) to be within the reach of MEG. The present bound already rules out a part of the parameter 
space. 

The SLIM model has a light sector (mass< 10 MeV) and a heavy sector (mass> mw)- The light sector consists of 
the scalar DM candidate, Si and at least one right-handed neutrino, N\. The heavy sector consists of the components 
of a scalar electroweak doublet: (i) a CP-odd neutral scalar (^2); (h) a CP-even neutral scalar (<5 2 ); (hi) a charged 
scalar (4>~). The light sector can show up in supernova explosion and the decay of light mesons. The present bounds 
from meson decay as well as supernova explosion are too weak to constrain the model. However, the on-going KLOE 
experiment [l3| and future NA62 experiment (l5j can test the model. In case of future observation of supernova 
explosion, invaluable information on this bound can be derived. Considering the fact that there is a lower bound on 
the coupling of the light sector to active neutrinos and an upper bound on the masses of the light sector, this model 
is falsifiable by low energy experiments with enough luminosity. 

The heavy sector can be produced at the LHC via the electroweak interactions. The present lower bounds are 
m - > 80 GeV [13 and ms 2 ,m c f >2 > 90 GeV For relatively light new particles (m^- < 130 GeV), LHC with 

30 fb _1 of integrated luminosity and at 14 TeV center of mass energy can make discovery via the pp —> 4> + 4>~ — > 
+ missing energy mode. In principle, by studying $T — > l~ + missing energy, the coupling \gi a \ 2 (see Eq. ([!])) 
or for the case that toat 2 <C m^- , \gi a \ 2 + \g2a\ 2 can be extracted. These are the same couplings that determine the 
neutrino mass matrix and the pattern of LFV rare decays. However, as discussed in to make this possible the 
large uncertainty induced by the uncertainties in the parton distribution functions should be reduced or the energy 
of center of mass should be increased to enhance the signal to background ratio. 

The other model that we discussed (AMEND) does not necessarily contain a low energy sector. Within this model, 
a scalar electroweak triplet (A) and two fermionic doublets with opposite hypercharges exist. A contains a doubly 
charged component, A ++ . This particle can lead to signals consisting of a pair of same sign charged leptons plus 
missing energy. The background to this signal from the SM is not very high so the discovery chance of the signal 
should be high provided that A ++ is not too heavy. 

These two models have quite different predictions for direct DM searches. Within the SLIM model, we do not 
expect a signal in the future DM searches. However, within AMEND, we generally expect a signal in the future DM 
searches. Only in some very specific parts of the parameter space such as when the DM mass approaches half the 
Higgs mass, the DM nucleon cross section is suppressed. The present bounds from direct DM searches already rule 
out a part of the parameter space. 
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